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Quarks produced in the early stage of noncentral heavy-ion collisions could develop a global spin 
polarization along the opposite direction of the reaction plane due to the spin-orbital coupling via 
interaction in a medium that has finite longitudinal flow shear along the direction of the impact 
parameter. We study how such polarization evolves via multiple scattering in a viscous quark-gluon 
plasma with an initial laminar flow. The final polarization is found to be sensitive to the viscosity 
and the initial shear of local longitudinal flow. 



PACS numbers: 25.75.-q, 13.88. +e, 12.38.Mh 

I. INTRODUCTION 

The observed jet quenching and collective phenom- 
ena in high-energy heavy-ion collisions at the Relativis- 
tic Heavy Ion Collider (RHIC) provide strong evidence 
of the formation of strongly coupled quark-gluon plasma 
(QGP) [1, 2]: The strong quenching of high transverse 
momentum jets is understood to be caused by parton 
energy loss induced by multiple collisions of the lead- 
ing parton with color charges in the thermal medium [3- 
8]; the observed collective flow in the final bulk hadron 
spectra indicates a hydrodynamic behavior of the initial 
dense matter as an almost perfect fluid with a very small 
shear viscosity [9, 10], r)/s < 0.5. The large jet transport 
parameter from the observed strong jet quenching and 
small shear viscosity inferred from the collective flow can 
be connected to each other through a transport process 
in a strongly coupled system [11]. They both describe 
the ability of the medium partons to transfer momentum 
via strong interaction in QCD and maintain local equi- 
librium. Globally, such transport processes help to dissi- 
pate variations of flow velocities and thus will reduce the 
anisotropic flow, which is driven by the initial geometric 
anisotropy [9, 10]. In this paper, we discuss the possi- 
bility of global quark spin polarization caused by such 
transport processes in noncentral high-energy heavy-ion 
collisions. 

It was first proposed by Liang and Wang [12] that 
global quark polarization could occur in the QGP formed 
in a noncentral heavy-ion collision. They argued that at 
a finite impact parameter, the initial partons produced 
in the collision can develop a longitudinal fluid shear dis- 
tribution representing local relative orbital angular mo- 
mentum (OAM) in the same direction as the global OAM 
of the noncentral nucleus-nucleus collisions. Since inter- 
action via one-gluon exchange in QCD contains a spin- 
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orbital coupling, the OAM could cause a global spin- 
polarization of quarks and antiquarks in the direction 
parallel to the OAM. Such a global (anti)quark polariza- 
tion should have many observable consequences such as 
global hyperon polarization [12, 13], vector meson spin 
alignment [12, 14], and the emission of circularly polar- 
ized photons [15]. Predictions have been made [12, 14-17] 
for these measurable quantities as functions of the global 
quark polarization P q . Experimental measurements of 
the A hyperon polarization with respect to the reaction 
plane at RHIC [18-25] place a limit |P A A | < 0.02 [19, 24]. 
Such a limit puts a stringent test on both the initial shear 
of longitudinal flow in noncentral heavy- ion collisions [17] 
as well as the time evolution of the polarization through 
transport processes. 

The estimates of the global quark polarization in 
Ref. [12] and in subsequent studies [16, 17, 26, 27] were all 
obtained by considering the polarization process for a sin- 
gle scattering between quarks and thermal partons. How- 
ever, one should consider the effect of the multiple scat- 
tering and expect that the quarks will be progressively 
polarized through multiple scattering. Furthermore, with 
the minimum values of shear viscosity rj/s > 1/47T in 
QGP imposed by the quantum limit, the local momen- 
tum shear, dp z /dx, of the fluid, that is, the local OAM 
of interacting parton pairs, will decay with time. This 
will lead to a nontrivial time evolution of quark polar- 
ization P depending on the shear viscosity of the QGP 
matter and the final state observed global polarization 
could serve as a viscometer of QGP. In this paper, we 
focus on these two issues with a simple and yet interest- 
ing hydrodynamic evolution of a relativistic laminar flow 
between two frictionless impenetrable walls. 

The rest of the paper is organized as follows. In Sec. II, 
we extend the calculation in Ref. [12] to the case of scat- 
tering of an initially polarized quark in a static potential 
model. In Sec. Ill, we study the relativistic laminar flow 
and compute the decay of the longitudinal momentum 
gradient. The results of Sec. II and Sec. Ill are applied 
to Sec. IV to study the time evolution of the quark po- 
larization. 
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II. POLARIZATION OF INITIALLY 
POLARIZED QUARKS 




I 



y' 




Reaction plane 


1 

















FIG. 1. (Color online) Illustration of noncentral collisions 
with impact parameter b of two heavy nuclei with radii Ra- 
The global angular momentum of the produced matter is 
along — y, opposite to the reaction plane. 



We consider two colliding nuclei with the beam pro- 
jectile moving in the direction of z and the impact pa- 
rameter b defined as the transverse distance of the pro- 
jectile from the target nucleus along the x direction as 
illustrated in the upper panel of Fig. 1. The direction 
y defines the reaction plane, y oc z x x. The initial 
OAM of these two colliding nuclei is along the direction 
opposite to the reaction plane and could be very large. 
Given 1 fm < b < 10 fm, the initial OAM L ~ Aby/s/2 
is roughly 10 5 < L < 10 6 for Au-Au collisions at RHIC 
energy Vs = 200 GeV and 3 x 10 6 < L < 3 x 10 7 
for Pb-Pb collision at Large Hadron Collider (LHC) en- 
ergy t/s = 5.5 TeV. Because of the unequal local num- 
ber density of the participant projectile and target nu- 
cleons at various transverse positions, some fraction of 
this large OAM could be transferred into the produced 
QGP matter in the overlapping region. Such global angu- 
lar momentum, however, would never lead to a collective 
rotation of the system since there is no strong binding 
or attractive interaction in the partonic interaction at 
high energy. Instead, it could be manifested in the fi- 



nite transverse (along x) gradient of the average longi- 
tudinal momentum p z per produced parton due to the 
partonic interaction at high energy (see the lower panel 
of Fig. 1). Given the range of interaction Ax, two col- 
liding partons will have relative longitudinal momentum 
Ap z = Axdp z /dx with relative OAM L y ~ — AxAp z . 
This relative OAM will lead to global quark polariza- 
tion along — y through the spin-orbital coupling in QCD. 
This is essentially the argument that was first proposed 
in Ref. [12]. It was found that the quark polarization via 
a single scattering with given relative momentum p reads 



P = 



Act 



cr-C + <7i 



TTfip 



2E(E + m) ' 



(2.1) 



where a s , s =t,| is the cross section of final quark with 
spin s along y, m is the mass of interacting quark, and 
(i is the Debye screening mass of longitudinal gluon, 
p? = g 2 (N c + N f /2)T 2 /3. The initial relative momen- 
tum p can be estimated as p ~ Axdp z /dx with Ax ~ /i -1 
being the characteristic range of interaction. Then pj \x 
is nothing but the relative orbital angular momentum 
between the scattering quarks, L y ~ —p//J-- In the non- 
relativistic limit for massive quarks, P is proportional 
to the spin-orbital coupling energy P oc E^s/p where 
E LS = (L ■ S){dV /dr)/rm 2 and (dV /dr)/r - ^ 3 with 
typical interaction range r ~ 1/yU. 

The estimates in Refs. [12, 16], and [17] were based on 
the assumption that the initial quarks are not polarized. 
In order to discuss the time evolution of the quark po- 
larization via multiple scattering, one must calculate the 
quark-quark cross section of initially polarized quarks. 
Let the fraction of initial quarks of spin A^/2 along y be 
R\ i = (l + AjPi)/2 with Pi being the initial polarization. 
The identity i?+ + i?_ = 1 must hold. Consider a quark 
with initial relative four- momentum p^ = (E, p) and spin 
Xi /2 scattering with a virtual gluon and resulting in final 
spin A//2; the cross section with fixed impact parameter 
xt is 
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(2.2) 



where Ct = 2/9 is the color factor associated with the 
target, (k-r) is the transverse momentum transfer 
from the virtual gluon to quark, and Pg^f.) ^ s ^ e nna l four- 
momentum of quark, P^r k \ = + [0, qr(kr)]. We use 
the screened static potential model to calculate M.x f \i 
in which = (A , 0) with A°(q T ) = g/(q T + fi 2 ) [3]. 

For small angle scattering (which is justified when the 
relative longitudinal momentum p is large), qT,kx ~ 
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From Eq. (2.2) and Eq. (2.3), it is evident that the po- 
larization will not change if one averages the cross section 
over all the possible directions of the parton impact pa- 
rameter xt- However, in noncentral heavy- ion collisions, 
the local relative OAM L y provides a preferred average 
reaction plane for parton collisions. This will lead to a 
global quark polarization opposite to the reaction plane 
of nucleus- nucleus collisions. This conclusion should not 
depend on our pcrturbative treatment of parton scatter- 
ing as far as the effective interaction is mediated by the 
vector coupling in QCD. Therefore, we average over the 
upper half-xy-planc with x > 0, that is, average over the 
relative angle between parton x-r and the nuclear impact 
parameter b from — tt/2 to tt/2 and over xt- To do this, 
we use the identity 



d 2 x T e i(kT - qT) - Xr 



2iri5{k y — q y ) 



i.c ••) k x - q x + i0+ 

Then the total unpolarized cross section reads, 



(2.4) 
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where K(x) is the complete elliptic integral of the first 
kind. The final global quark polarization is then 



Pf=Pi 



(i - p 2 y m 

2E(E + m) - Pi-Kfip' 



(2.7) 



It is also useful to get the transverse momentum de- 
pendence of the quark polarization. From Eq. (2.5) and 



Eq. (2.6), we read out the differential cross sections, 
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The transverse-momentum-dependent polarization 
(TMDP) defined as P f {q T ) = {dAa / dq T ) / {da / dq T ) now 
reads, 



Pfi.tr) 



nE{E + m)P l - p^ 'q 2 + \x 2 K{q T / V# + » 2 ) 



irE{E + m) - P % p^Jq\ + li 2 K{q 1 
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Some discussions are in order. (1) If the initial quark is 
unpolarized, Pi = 0, we recover the result of Ref. [12]. (2) 
Because the denominator is always positive in the right- 
hand-side (RHS) of Eq. (2.7) for high relative longitudi- 
nal momentum (i.e., when small angle approximation is 
applicable) , we always have Pf ^ Pi. Therefore scattered 
quarks always prefer to be polarized along — y direction. 

(3) The scattering matrix elements ^\ s \i with spin flip- 
ping (Xf = —Xi) are zero according to Eq. (2.3), so there 
is no flipping of quark's spin via the scattering under this 
small angle approximation. The polarization in the final 
state is caused by the larger cross section of quarks with 
spin up relative to quarks with spin down. This will lead 
to the conclusion that if the initial quark is completely 
polarized, Pi = ±1, we must have Pf = Pi. This is 
indeed the case expressed in Eq. (2.7) when Pj = ±1. 

(4) The quark polarization has a remarkable transverse 
momentum dependence, as shown in Eq. (2.10). Fig- 
ure 2 shows the typical behavior of TMDP as a function 
of the transverse momentum with given p = 10/x. The 
polarization grows with the transverse momentum due 
to quark-quark scattering. In principle, the A-hyperon 
polarization should have similar transverse momentum 
dependence, although as we mentioned in Sec. I it is not 
trivial to construct a correspondence between quark po- 
larization and hadron polarization. 



III. RELATIVISTIC LAMINAR FLOW 

Before discussing how the quark polarization evolves 
in a viscous QGP due to multiple scattering, we have to 
know how the QGP itself evolves through either trans- 
port model or viscous hydrodynamical model [9, 10, 28- 
32], Moreover, we also have to know the initial profile of 
the longitudinal flow field. In the discussion in Sec. II, 
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FIG. 2. The TMDP as a function of transverse momentum 
in unit of fj,. The initial relative longitudinal momentum is 
chosen to be p = 10/x. 
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FIG. 3. (Color online) Illustration of the velocity profiles of 
the relativistic laminar flow. 



we simply followed Ref. [12], and assumed that noth- 
ing depends on the longitudinal position in the system. 
In such a case, the finite angular momentum must lead 
to a velocity profile depicted in Fig. 3 (see Ref. [33] for 
a discussion of possible consequences of such a profile). 
On the other hand, another extreme is to assume that 
dvz/dx = everywhere, but the angular momentum is 
carried by the matter distribution in the reaction plane; 
see Ref. [34] for illustration. 

To study the effect of viscosity on the decay of the 
local angular momentum, we consider a simple laminar 
flow without driving force between two frictionless (free- 
slip flow) impenetrable walls. We assume the walls arc 
infinitely large and separated by a distance 2h. To make 
dimensions relevant for a heavy-ion collision, we set h = 5 
fm. Such a scenario might be far from the real longitu- 
dinal flow profile in high-energy heavy-ion collisions, but 
it will be very illustrative for our study here. We further 
assume that the flow profile has no longitudinal variation 
and the system has a reflection symmetry respect to the 
j/z-planc. We study two cases: One with no expansion, 



and another with boost-invariant expansion in y direc- 
tion, that is, with flow profile v y = y/t. In both cases, 
the flow four-velocity in the reaction plane has the gen- 
eral form = (7, jv x , 0, jv z ) with 7=1/ yl — u| — u| 
and v XjZ (t, x) being the x and z components of the three- 
velocity. 

As is well-known, the relativistic Navier-Stokes hy- 
drodynamics is unstable and provides a possibility for 
acausal signal velocities [35]. Therefore, we use the 
second-order theory by Israel and Stewart [36] instead. 
Although hydrodynamics has been widely used to model 
the heavy-ion collisions, as far as we know, there is no 
literature discussing the relativistic laminar flow. 

If there are no conserved charges, the hydrodynami- 
cal equations of motion are given by the conservation of 
energy and momentum 

9 M T^ = 0, (3.1) 

where T^" = (e + Q)u"u v - Qg^ + ir^ is the energy- 
momentum tensor, e is the energy density, is the pres- 
sure 1 , and 7r Miy is the shear stress tensor. To close the 
set of differential equations, one also needs to specify an 
equation of state (EOS) e = e(0). For simplicity, we use 
the ideal gas EOS e = 36. 

In its simplest form, Israel-Stewart hydrodynamics 
means that instead of being directly proportional to the 
velocity gradients, the shear stress tensor is a dynamical 
variable, which relaxes toward the Navier-Stokes value 
on its relaxation time r v : 

Dtt^ = ~ (n^ - 2r)V^u u} ) - 2n K ^u l,) Du K {3.2) 

where D = u x d x , A^ = (A^ + A^)/2, A^ = 
[A^A^ - ±A' iV A a p]A a P , V = - u^u v d v , A» v = 
— u^u v and r\ is the shear viscosity coefficient. The 
last term is required to keep the shear stress tensor or- 
thogonal to the flow velocity in all circumstances. This 
is the so-called truncated Israel- Stewart equation. Al- 
though there are more terms in a complete Israel- Stewart 
equation, for our purpose here, the truncated one is suf- 
ficient. 

The relaxation time is given by [36] 

TV = 2^2, (3.3) 

which is dependent on the shear viscosity and another 
coefficient fa. For massless Boltzmann particles, the ki- 
netic theory gives [36] 



If there is no phase transition, it is expected that fa for 
Fcrmion and Boson gases have only minor modification 



1 Note that since we used P to denote polarization, to avoid con- 
fusion, we do not use it to denote pressure. 
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from /?2 for Boltzmann gas at high temperature [37-42]. 
Taking temperature T ~ 350 MeV, the relaxation time is 
around r ff ~ 0.27—1.35 fm if using -q/s = l/(An) — 5/ (47r), 
where s is the entropy density, and for free gluon gas it 
is 



27r 2 r 3 



(3.5) 



with the degeneracy factor v g — 2(N^ — 1). 

Since the system has reflection symmetry with respect 
to the yz-plane, and there are no particle, momentum, or 
heat flow through the hard walls, the system obeys the 
following boundary conditions 



««(*,0) =v x (t,0) =v x (t,±h) =0, 
dv z {t,±h)/dx = 0. 



(3.6) 



As the initial state we choose uniform initial temperature 
of 355 MeV (corresponding to RHIC initial temperature) , 
no flow in ic-direction, and a simple sine-type longitudinal 
flow velocity profile 



v x {t ,x) = 0, 
v z (t , x) = vq sin (nx/2h), 



(3.7) 



where vq is the magnitude of the initial velocity at the 
two boundaries. In the following numerical calculation 
we consider two cases: vq = 0.7 and 0.9. In the ex- 
panding case, we use the initial time To = 1 fm. Since 
the shear stress tensor is a dynamical variable in Israel- 
Stewart hydrodynamics, we need its initial value too. A 
natural choice is the Navier-Stokes value, but its exact 
evaluation is difficult. It contains the time derivative of 
the flow velocity, which is unknown before the hydro- 
dynamic equation is solved. To avoid this problem, we 
initialize the shear stress, not to its exact Navier- Stokes, 
but to a "static Navier-Stokes" value; that it, we ignore 
all the time derivatives in the Navier-Stokes definition of 
the shear stress tensor and calculate the value based on 
spatial derivatives only. In practice this means that some 
components of the tensor are slightly larger and some 
slightly smaller than their exact Navier-Stokes values. 

In Fig. 4, we depict the time evolution of the gradient 
of the longitudinal momentum per particle averaged over 
x G [0, h], 



dpz 
dx 



h dxJ°(x)±^ 
XJ [X) dx J°(x) , 



dxJ°(x) 



(3.8) 



where J° = 7/3 is the proper particle number density. As 
expected, the shear viscosity dissipates the average gra- 
dient of the longitudinal momentum, especially for larger 
values of shear viscosity. The transverse expansion accel- 
erate this degradation, since strong transverse expansion 
means larger shear (shear tensor). 

In the case of transverse expansion and large viscosity, 
there appear to be a "shoulder" in the time evolution of 
the longitudinal momentum gradient (dp z /dx) as shown 
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FIG. 4. (Color online) Evolution of the average gradient of 
the longitudinal momentum per particle, dp z /dx, at different 
shear viscosities. Upper panel: the system has no transverse 
expansion. Lower panel: the system has Bjorken expansion 
in the y direction. 



in the lower panel of Fig. 4, where the gradient drops very 
fast initially and then slows down for a while before it de- 
creases again. The temporary slow down is caused by the 
oscillatory behaviors of the induced transverse flow in in- 
direction, and the particle number density J , which is 
used as a weight in the calculation of the average longitu- 
dinal momentum gradient in Eq. (3.8). The oscillations 
are an artifact of the fixed wall boundary conditions in 
our simple scenario. When there is no transverse expan- 
sion, the degradation is slower and there is no shoulder 
because of the smaller shear (in shear tensor). 

In Fig. 5, we show the profiles of velocity v z at different 
time with viscosity t]/s = 5/4-7T with (lower panel) and 
without (upper panel) transverse expansion. One of the 
functions of the shear viscosity is to transform the kinetic 
energy of the fluid to internal energy, hence damping the 
fluid shear (as shown in Fig. 5) and heating up the fluid. 
This can be explicitly seen in the upper panel of Fig. 6, 
where the temperature evolution is shown for the non- 
expanding system. The transverse Bjorken expansion in 
our problem, however, will dilute the system and cool 
the system down, overcoming the slight heating-up by 
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FIG. 5. (Color online) The profile of longitudinal velocity v z 
at different times with rj/s — 5/4-7T. 
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following evolution equation for the polarization, 



the shear viscosity, as shown in the lower panel of Fig. 6. 
The transverse expansion will also accelerate the degra- 
dation of the longitudinal velocity as shown in the lower 
panel of Fig. 5 as compared to the upper panel for the 
case of no transverse expansion. 



IV. EVOLUTION OF THE GLOBAL QUARK 
POLARIZATION 



With the model of time evolution of the longitudinal 
momentum gradient of the medium partons we can now 
study the time evolution of the quark polarization when 
it is progressively polarized due to multiple scattering. 

According to Eq. (2.7), the change of polarization 
caused by one scattering is 



dP _ AP 1 (1 - P 2 )TTfl P 

~dt = ~ ' ~T q ~2E(E + m)-Pnpp' 



(4.2) 



where r q is the mean-free-path of quark which is related 
to the transport cross section a t r of the interacting par- 
tons though T q ~ l/(pa t r), where p = is g ((3)T 3 /n 2 is 
the density of medium gluons, assuming gluons are the 
dominant degrees of freedom in the medium. The shear 
viscosity for a thermal ensemble of gluons is roughly [43] 



1 



n 



T-pr q . 



(4.3) 



We have then the final rate equation for the time evolu- 
tion of the quark polarization, 



dP 
~dt 



4Tps (l-P 2 )irpp 



9s r)2E{E + m) - Pnpp' 



(4.4) 



AP = Pt — Pi = - 



(1 - P 2 )irpp 



2E(E 



P l TXpp 



(4.1) 



For convenience we denote P = Pi. Then we get the 



From Eq. (4.4), the rate dP/dt is inversely propor- 
tional to the viscosity. This is evidently shown in the 
upper panel of Fig. 7, in which the evolutions of the 
quark polarizations arc shown for the initial polarizations 
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FIG. 7. (Color online) Evolution of the average polarization 
P — Aa/a with initial polarization P(to) = with different 
values of viscosities without (upper panel) and with (lower 
panel) transverse expansion. 



P(0) = and for a system without transverse expansion. 
With transverse expansion, the mean-free-path increases 
more rapidly with time and therefore slows down the po- 
larization rate. The transverse expansion also accelerates 
the degradation of the longitudinal momentum gradient, 
reducing the polarization in each scattering. Both effects 
slow down the time evolution of the polarization in an 
expanding system as shown by comparison between the 
upper and lower panels of Fig. 7. 

Because of the reheating by viscous interaction, the 
initial cooling of the system due to transverse expansion 
is significantly slower for a larger value of shear viscos- 
ity, as shown in Fig. 6. This speeds up the polarization 
according to Eq. (4.4). However, a larger shear viscosity 
also slows down the polarization because the polariza- 
tion rate is inversely proportional to the shear viscos- 
ity. During the early stage of evolution, the second effect 
dominates, leading to a slower polarization process with 
a larger value of shear viscosity. At a later time, effect 
of reheating becomes more dominant and a larger shear 
viscosity leads to a faster polarization process. 

The polarization is also sensitive to the initial condi- 



tion of the longitudinal flow shear. In our simple laminar 
flow model, the initial longitudinal flow shear is propor- 
tional to the value of vq. The nonlinear dependence of the 
polarization rate on the relative momentum p in Eq. (4.4) 
determines the nontrivial dependence of the polarization 
on the values of vq as shown in Fig. 7. 

Note that the polarization rate we used are derived 
with the approximation of small angle scattering which 
is only valid when the longitudinal momentum gradient is 
large. For large shear viscosity r\j s and at late time, the 
longitudinal momentum gradient can become too small. 
One can no longer use the rate equation derived here. 
However, one can assume that the polarization process 
will stop at this point when there is not significant local 
orbital angular momentum. 



V. SUMMARY 

In conclusion, we have calculated the polarization cross 
section for quarks with initial polarization within the 
frame of perturbative QCD, which we use to study the 
time evolution of the quark polarization via multiple scat- 
tering in a medium with nonvanishing local orbital angu- 
lar momentum. We considered the simple case of laminar 
flow as governed by viscous hydrodynamics with given 
shear viscosity rj/s and a simple illustrative initial con- 
dition. Such a simple hydrodynamic model provides the 
dynamic evolution of the longitudinal flow shear as the 
polarization mechanism for quarks via parton scattering. 
Because the values of the shear viscosity influence the 
degradation of the longitudinal flow shear with time and 
the cooling of the system, it also determines the time 
evolution of the quark polarization. Since the polariza- 
tion rate is inversely proportional to the shear viscosity 
and depends nonlincarly on the average longitudinal mo- 
mentum shear, the final quark polarization is found to 
be sensitive to the shear viscosity but has a nontrivial 
dependence. In this sense, one can use the final state 
polarization as a possible viscometer of the QGP. 

For more realistic studies, one should employ a full 
scale 3+1D viscous hydrodynamics [32] with initial con- 
ditions from Monte Carlo models such as HIJING [44]. 
The initial parton production from this kind of model 
has approximate Bjorkcn scaling which will give rise to 
very small initial local longitudinal flow shear [17] except 
at very large rapidity region. Such small initial local 
longitudinal flow shear comes from the violation of the 
Bjorkcn scaling which one can use as the initial condi- 
tion. Furthermore, one should also extend the current 
calculation of the quark polarization beyond the small 
angle approximation. 
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